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Abstract 

A method for detecting the transition Vr in long-basehne accelerator 

experiments, that consists in comparing the far-to-near ratios of the spectra 
of quasielastic CC events generated by high- and low-energy beams of muon 
neutrinos, is proposed. The test may be accessible to big water Cherenkov 
detectors and iron-scintillator calorimeters, and is limited by statistics rather 
than systematics. 
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The data of Super-Kamiokande favor the transition z/^ — * Vj. as the source 
of the deficit of muon neutrinos from the atmosphere. However, this still has to be 
verified by directly observing v^- appearance in accelerator long-baseline experiments. 
The options discussed thus far, all involving fine instrumentation on a large scale, 
are to detect the secondary r by range in emulsion 0, by Cherenkov light of the 
T Q, or by the transverse momentum carried away by the decay neutrino (s) 0]. 
By contrast, in this paper we wish to formulate a r signature that is solely based 
on the energy spectra of CC events, and therefore should be accessible to water 
Cherenkov detectors and to relatively coarse calorimeters with muon spectrometry. 
We assume that the experiment includes a near detector of the same structure as the 
far detector, irradiated by the same neutrino beam but over a short baseline that 
rules out any significant effects of neutrino oscillations 0. Thereby, the systematic 
uncertainties in comparing the interactions of primary and oscillated neutrinos are 
largely eliminated. Our aim is to distinguish the muonic decays of r leptons against 
the background of i/^-induced CC events. In order to minimize the effects of 
disappearance, the data collected with a harder r-producing beam are compared 
with those for a softer reference beam in which r production is suppressed by the 
threshold effect. The analysis is restricted to quasielastics (QE), that is, to neutrino 
events featuring a muon and small hadronic energy. 

As soon as the first maximum of the oscillation lies below the mean energy of 
muon neutrinos in the beam, or Am?L/ [EJ) < 1.24 eV^km/GeV, much of the signal 
from QE production and muonic decay of the r is at relatively low values of visible 
energy E. That is because the tau neutrinos arising from are softer on 

average than muon neutrinos, the threshold effect is relatively mild for quasielastics, 
and a large fraction of incident energy is taken away by the two neutrinos from 

li~vv. Let /(-E) be the spectrum of QE events observed in the far detector, 
n(Ei) — the spectrum of similar events in the near detector that has been extrapolated 
and normalized to the location of the far detector, and R[E) — the ratio of the two: 
R{Ei) = f (E) / n{E) . In the case of z/^ disappearance through the transitions z/g 
or ^ where is the hypothesized sterile neutrino, the ratio R for the harder 
beam should be identically equal to that for the softer beam: i?hard(-E) = Rsoft{E). 
However, in the case of z/^ — > Ur this equation is violated by the process of r 
production and muonic decay, that predominantly occurs in the harder beam and 
shows up as a low-E enhancement of the corresponding "far" spectrum /hard(-E). 
This causes the ratio -Rhard to exceed i?soft towards low values of visible energy 
E. The latter effect, that may provide a specific signature of appearance, is 
investigated in this paper. 

In the simulations reported below, the harder (or r-producing) and softer 
(or reference) beams of muon neutrinos are respectively assigned as the high- and 
low-energy beams from the Main Injector (MI) at Fermilab, as foreseen by the 
NuMI program ||^. The mean z/^ energies in these beams are close to 12 and 5 
GeV, respectively. Again as in the NuMI program, a baseline of 730 km is assumed 
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throughout. The systematic uncertaities on the near spectra n{E) are neglected^. 
The two detector types considered are an iron-scintillator calorimeter and a water 
Cherenkov detector. 

Charged- current interactions of the z/^ and Ur are generated using the NEU- 
GEN package that is based on the Soudan-2 Monte Carlo 0. This generator takes 
full account of exclusive channels like quasielastics proper and excitation of baryon 
resonances, that are important for our analysis of CC events with small hadronic 
energy. 

The iron-scintillator calorimeter is assumed to be the MINOS detector 
that is in construction stage. The detector response is not simulated in full de- 
tail; instead, the resolution in muon energy is approximated as 6E^ = 0.11 x E^^^'^ 

and in energy transfer to hadrons — as Su (GeV) = 0.55 x 0. Quasielastic 

events are selected as those with E^^ > 800 MeV and v < 1 GeV, where E^ and v 
are the smeared values of muon energy and of energy transfer to hadrons, respec- 
tivelyQ. Given the characteristic topology of such events in the detector (a single 
track traversing more than three nuclear interaction lengths in iron plus a few scin- 
tillator hits near the primary vertex), we assume that they are reconstructed with 
100% efficiency and that the background from pion punchthrough is insignificant. 
The visible energy of a CC event, E, is estimated in terms of smeared quantities: 
E = E^ + p. 

We also analyze the performance of a water Cherenkov detector, in which v^- 
induced quasielastics have a characteristic signature of muonlike single-ring events 
[0. Analyzing the multi-GeV exiting muons, particularly in the smaller near detec- 



tor, may require a device of the AQUA- RICH type |10] or an instrumented muon 
absorber downstream of the water tank. In our simulation for the water Cherenkov 
detector, QE events are selected as those featuring a single muon with momen- 
tum above 200 MeV, no additional charged secondaries with momenta above the 
Cherenkov threshold in water, and no vr'' mesons in the final state. For the quasielas- 
tic reaction u^n fJ'~P, this implies an energy transfer to the nucleon of less than 
0.47 GeV. (Some 70% of thus selected events are due to quasielastics proper.) Muon 
energy is then smeared according to 6E^ = 0.05 x E^^'^'^, and neutrino energy is es- 
timated assuming quasielastic scattering on a free neutron: 

^_ mNE^-ml/2 

niN - Ef,+Pf, cos Of, ' 



^The systematic uncertainty on n[E), that arises from correcting the spectrum in the near 
detector for the neutrino source not being pointHke, is analyzed in |^. For the statistics consid- 
ered in this paper, the overall uncertainty on the effect is still dominated by statistical errors. 
Furthermore, we may expect that systematic uncertainties on n{E) for the high- and low-energy 
beams from the Main Injector are correlated and, therefore, will partially cancel in the difference 
of far-to-near ratios for these two beams. 

^ These selections should be viewed as illustrative. The actual selections will be based on a 
detailed simulation of detector response to CC events with small hadronic energy. 
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Here, mjy and are the neutron and muon masses, and and are the muon 
momentum and angle relative to incident neutrino, respectively. 

Shown in Fig. |^ are the oscillation-free near spectra of selected QE events, 
n{E), for the two beams and two detectors considered. In the absence of oscillations, 
10 (50) kton-year exposures of the iron-scintillator (water Cherenkov) detector in 
the softer and harder beams will yield some 1200 (1600) and 4700 (4900) z/^-induced 
QE events, respectively^. Assuming either — > 1/^ or driven by Am^ = 

0.01 eV^ and maximal mixing of sin^2^ = 1, the far-to-near ratios R{E) for either 
beam and either detector are illustrated in Fig. ^ That the ratios R{E) for the 
transitions — > z/^ and z/^ — > Vs diverge towards low values of E is evident for the 
harder beam in which r production is not suppressed. Again considering 
and z/^ Us with maximal mixings, in Fig. ^ we plot the difference 

AR{E) = -Rhard — -Rsoft 

for various values of Am?. Indeed, at visible energies below some 4 GeV AR{E) 
deviates from zero for the transition z/^ ^ z/^, while staying close to zero for ^ Vs- 
This deviation may be viewed as a signature of Vr appearance. The naive expectation 
for —y z/g, AR{E) = 0, is violated by the smearing of neutrino energy, and 
therefore is better fulfilled for the water Cherenkov detector. 

By the time the proposed test can be implemented, the actual value of Am^ 
will probably be estimated to some 10% by analyzing the z/^ disappearance in the 
MI low-energy beam |p. Given the value of Am^, a consistent approach would 
be to fit AR{E) to the predicted shape in order to estimate the mixing between 
the muon and tau neutrinos. A cruder measure of the effect is provided by the 
integral S = J AR{E)dE, which we estimate between E = 1 and 3 GeV for the 
iron-scintillator detector, and between E = 0.5 and 3 GeV for the water Cherenkov 
detector. Allowing for either z/^ and z/^ Ug with maximal mixing, the 

respective integrals 5'(z/^ Vr) and S{u^ u^) are plotted in Fig. ^ as functions 
of Am?. 

Since the reference beam produces many more low-energy events than the 
harder beam, the statistical error on the integral S is largely determined by the 
statistics accumulated with the latter. Therefore, we fix the exposure of the iron- 
scintillator (water Cherenkov) detector in the reference beam at 10 (50) kton-years, 
and assume a similar or bigger exposure in the harder beam: 10(50), 20(100), or 
40(200) kton-years. The respective statistical uncertainties on ^(z/^ Vr) are 
illustrated by successive error corridors in Fig. ^ And finally, dividing the difference 
between ^(z/^ Vr) and ^(z/^ z/g) by the statistical error on S{u^ u^), we 

^In all numerical estimates, we do not take into account the discussed upgrade of the proton 
driver at Fermilab that may result in a substantial increase of neutrino flux from the Main 
Injector. 

^In the Figures, statistical fluctuations are suppressed for the data points themselves, but the 
error bars are for the statistics as indicated. 
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estimate the statistical significance of tfie enliancement that is also depicted in Fig. 
I- 

We estimate that at a level of 3a, the 10(50), 20(100), and 40(200) kton- 
year exposures of the iron-scintillator (water Cherenkov) detector in the MI high- 
energy beam will allow to probe Vr appearance down to the Am^ values of some 
0.008 (0.005), 0.006 (0.004), and 0.005 (0.003) eV^, respectively. Thus in the NuMI 
program with the 5.4-kton MINOS detector and with the existing Proton Booster 
0, the proposed test may be sensitive to Am^ values in the Kamiokande-allowed 
region |T2[, but not below some 5 x 10~^ eV^ as suggested by the more recent results 



of Super-Kamiokande |jT3[- On the other hand, a big (~ 100 kton) water Cherenkov 
detector would allow to probe the transition over a large portion of the 

Am^ region favored by the analysis of atmospheric neutrinos in Super-Kamiokande[|. 

To conclude, we have proposed a test of appearance that consists in com- 
paring the far-to-near ratios of the spectra of quasielastic CC events generated by dif- 
ferent beams of muon neutrinos, and therefore may be accessible to water Cherenkov 
detectors and to calorimeters with muon spectrometry. The test is limited by statis- 
tics rather than systematics, and its significance crucially depends on the exposure 
in the harder beam in particular. 
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Figure 1: The oscillation-free "near" spectra of i/iz-induced quasielastic events, n{E), 
for the iron-scintillator (on the left) and water Chcrenkov (on the right) detectors 
irradiated by the low-energy and high-energy beams from the Main Injector. The 
assumed exposures are 10 and 50 kton-years, respectively. 
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Figure 2: Assuming the transitions z/^, and z/g (solid and open dots, 

respectively) driven by Am^ = 10^^ eV^ and sin^ 26 = 1, the far-to-near ratio 
R{E) = f{E)/n{E) for quasielastic events produced by the MI high- and low- 
energy beams (top and bottom panels, respectively) in the iron-scintillator and 
water Cherenkov detectors (left- and right-hand panels, respectively). The error bars 
on R{E) for the transition — > z/^ are the statistical uncertainties corresponding 
to 10 (50) kton-year exposures of the iron-scintillator (water Cherenkov) detector 
in either beam. 
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Figure 3: Assuming either z/^ — >■ Vr or — > z/g (solid and open dots, respectively) 
with maximal mixing, the difference AR{E) between the far-to-near ratios for the 
MI high- and low-energy beams. Simulated data for the iron-scintillator and water 
Cherenkov detectors are shown in the left- and right-hand panels, respectively. The 
top, middle, and bottom panels are for Am^ = 0.004, 0.007, and 0.010 eV^, respec- 
tively. Depicted by error bars are the statistical errors that correspond to 10 (50) 
kton-year exposures of the iron-scintillator (water Cherenkov) detector in either 
beam. 
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Figure 4: Assuming either z/^ — ^ z/^ or z/^ ^ z/g (solid and open dots, respectively) 
with maximal mixing, the integrated difference S (see text) as a function of Am^ for 
the iron-scintillator (top left) and water Cherenkov (top right) detectors. Shown 
by successive error corridors are the statistical errors on 5'(Am^) corresponding 
to 10(50), 20(100), and 40(200) kton-year exposures of the iron-scintillator (wa- 
ter Cherenkov) detector in the r-producing beam. The bottom panels show the 
statistical significance of the r signal for either case. 



10 



